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We study two body dipelar scattering with one dimension of confinement. We include the effects 
of confinement by expanding this degree of freedom in harmonic oscillator states. We then study 
the properties of the resulting multi-channel system. We study the adiabatic curves as a function 
of D/l, the ratio of the dipolar and confinement length scales. There is no dipolar barrier for this 
system when D/l < 0.34. We also study the WKB tunneling probability as a function of D/l and 
scattering energy. This can be used to estimate the character of the scattering. 

PACS numbers: 34.20.Cf,34.50.-s 
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I. INTRODUCTION 

There has been great progress in the production of 
ultracold polar molecules [l|, 0| and ultracold dipolar col- 
lisions have now been experimentally observed Q. The 

dipolar interaction is: Vdd = <P "'^"'^^a i where d is the 
induced dipole moment along the field axis (z). When 
this is included in the ultracold environment the scatter- 
ing properties are intriguing [^-Q- For example, dipolar 
coUisions are not Wigner suppressed at threshold; rather 
non-zero partial wave cross sections go to a Constant as 
energy goes to zero These ultracold molecular sys- 
tems present exciting avenues to study chemical reactions 
at ultracold temperatures with unprecedented control [1- 
[ip'|. There is a down side to chemical reactions, it may 
shorten trapping lifetimes and prevent the production of 
degenerate dipolar gases. To further control the collisions 
one could use confinement and produce lower dimensionai 
quantum gases. By placing the ultracold molecules in a 
ID optical lattice aligned with the polarization axis, one 
would remove the possibility of the molecules approach- 
ing each other along the attractive configuration of the 
interaction. In the strong trapping limit, this would make 
the interaction a purely repulsive cP/p^ potential. This 
would inhibit the particles from reaching the short range 
where they would chemically react. Beyond few body 
physics, ultracold 2D systems have been used to study 
the BKT phase transition [ll| and quasi-condensation 
[m . Adding dipoles to these systems could lead to exotic 
dipolar many body systems [13i] and quantum memories 

li. 



The production of 2D dipolar gases would be a valu- 
able addition to the study of ultracold gases. They will 
ofFer further control over collisions and enable production 
of exotic many body systems. In an effort to bolster their 
production, we study basic properties of the 2D dipolar 
system. We are most interested in understanding when 
there is a dipolar barrier and how resilient a barrier it 
is. Previous studies of ultracold 2D scattering have fo- 
cused on ultracold atoms jT^ and inelastic collisions [l^ . 
These theories were based on the assumption that the 



short range interaction is much smaller than the con- 
finement length. We will use the same assumption here, 
but we allow the dipolar length scale to be both large or 
small compared to the confinement length scale. Con- 
fined dipolar systems have been considered in ID dipolar 
[13], but to the knowledge of the author, no one has yet 
studied quasi 2D dipolar scattering. Here we study the 
scattering of two point dipoles under 1 dimension of con- 
finement. By considering channels formed by the eigen- 
states of the confinement Hamiltonian which are coupled 
to each other by the dipolar interaction, we are able to 
study a standard multi-channel scattering problem. With 
this model we address an important issue: when can the 
particles be prevented from reaching short range. This 
work shows when confinement can be used to effectively 
shield the molecules from short range chemical reaction 
processes 0]. In addition, we offer an estimate of char- 
acter of the scattering, i.e. threshold or semi-classical. 
This will be important to understand the interactions 
and the properties of the resulting many body system. 
Furthermore, we address whether or not this system be- 
haves similarly to a pure 2D system. This comparison 
allows one to predict the 2 body scattering properties of 
the system easily. The rest of the paper contains a review 
of the scattering equations, then a study of the adiabatic 
curves, and finally a study of the WKB tunneling proba- 
bility. This allows rough scattering characteristics to be 
determined about the 2D dipolar system. 



II. EQUATIONS OF MOTION 

The Schròdinger Equation for two dipoles with 1 di- 
mension of confinement is 
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Here we assume the particles are polarized along the 
confinement axis (z). The length scales of this system 
are the dipolar length D = fid? /hj^ and the confinement 
length l = ^yh/fiuj where fi is the reduced mass and oj 
is the frequency. The corresponding energy scales are 
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FIG. 1: (Color Online) (a) The adiabatic potential energy 
curves, Vo{p/l), are shown for D = 0.5 (solid black), D = 1 
(red circles), D = 2 (blue x), and D = 5 (green square) 
in harmonic oscillator units. (b) The height of the dipolar 
barrier as a function of D. The first excited threshold, 2D^ 
(solid black) and D^''^ (blue dashed) are plotted in dipolar 
energy units. 

Eu = fi^ /2fj.D'^ and hu. Additionally, there is the scat- 
tering energy E, and this has a length scale of k^^ where 

~ 2iiE/h?. The parameters of interest in this system 
will be the ratio of the length scales D = D/l and the 
scattering energy. By varying D, one will achieve the 
coUisional control. 

To study this system, we expand the wavefunction in 
harmonic oscillator states in the confinement direction, 
Xn(z), and partial waves in the azimuthal coordinate, 
girrup^ Using these expansions the total wavefunction is: 
V' = Em„ e™'^x„(z)^i^ [3, and this allows us to in- 
tegrate out z. Now using D to rescale Eq. ([T]), we obtain 
a multi-channel radiai Schròdinger equation describing 
quasi 2D dipolar scattering: 

f 1 £ m2 - 1/4 -2 E \ ^ , . 

-^Unn'{p)(l>7nn'{p) (2) 
mn' 

where p — p/D and 

{Xn\Vdd\Xn') = d^Unn'ip) 

Unn'ip) = J dzxliz) ^^^^ ^^^2)3/2 ^"'(^)- 

We set the lowest threshold to zero and consider only 
even or odd n because parity is conserved by the inter- 
action. The short range length scale is assumed to be 
much smaller than both l and D. This length scale is 
where the interaction deviates from the dipolar interac- 
tion, and is typically of order lOOao (^5nm). This length 
scale will most likely be the depolarization length where 
the molecules reorient relative to each other. Once the 
particles reach this length scale more complex interaction 
dynamics occur, such as chemical reactions and are far 
beyond the scope of this work. 

To study the character of the system, we first look at 
the adiabatic curves. These are obtained by diagonalizing 
the potentials in Eq. ([5]) as a function of p. We also 
include the standard diagonal correction from the kinetic 




FIG. 2: (Color Online) (a) The lowest adiabatic curve on a 
log-log scale for various values of D, which are: 0.5 (black 
triangle), 1 (red circles), 2 (blue x), 5 (green square), and 
10 (violet diamond). Compared with the pure 2D scattering 
potential (solid black). (b) The lowest diabatic curve is shown 
for _D=0.1 (blue -I-), 0.34 (black triangle), and 1 (red circles) 
as a function of p. Also shown is pure dipolar repulsion (solid 
black) and the absolute value of the attractive 2D centrifugai 
terni (solid blue). 



energy which is -^{0\(P /dp'^\0) where |0) is the lowest 
adiabatic eigenstate and depends on p. First we will look 
at the adiabatic potentials in confinement units. In Fig. 
[1] (a) the lowest m = adiabatic potential energy curves 
are shown for D — 0.5 (solid black), D — 1 (red circles), 
D = 2 (blue x), and D — 5 (green square). This figure 
shows that as D is increased the barrier becomes higher 
and moves out in p/l. In Fig. [1] (b) we have shown 
the height of the adiabatic barrier as a function of D 
(solid red) in dipolar units, which more readily translate 
to scattering physics. We have also plotted the energy 
of the first excited threshold, 2D^ (solid black) and D^^'^ 
(blue dashed) [l^. In this figure the barrier height goes 
to zero at non-zero D. This will be explored below. 

To relate these adiabatic potentials to the scattering 
character it is convenient to look at them in dipolar units; 
in these units the long range interaction has a simple 
form: (m^ — l/4)/2p^ + l/p^. The scattering on this po- 
tential is understood and has a universal form [T^ . Thus 
it will be important to see when the adiabatic curves de- 



viate from the pure 2D potential and this will offer an 
idea of when the scattering deviates from the pure 2D 
case. 

The adiabatic potentials in the dipole units are shown 
in Fig. [5] (a) for D equal to: 0.5 (black triangle), 1 (red 
circles), 2 (blue x), 5 (green square), and 10 (violet di- 
amond), and the pure 2D case is sohd black. For large 
values of D, there is a large dipelar barrier. The ki- 
netic energy correction makes a significant contribution 
for large D. Importantly the dipolar barrier closely mim- 
ics the pure 2D case. As D is decreased, the region where 
the potential deviates from the repulsive 1//?^ moves out 
in p. This deviation from the pure 2D case is very clear 
at _D = 0.5 (black). If D is further decreased, eventually 
there is no barrier below 0.34 (= De). 

In Fig. [5] (b) , we illustrate why there is no barrier be- 
low De, and it is essentially a signal channel phenomena. 
We plot the diabatic curves, Uoo{p), for different values 
of D^O.l (blue +), 0.25 (green dashed), 0.5 (black tri- 
angle), and 1 (red circles) as a function of p. The pure 
2D repulsion (l/p^) is shown as a solid black line and 
the absolute value of the attractive centrifugai term is 
show as solid blue (1/8/5^). At small values of p/l, the 
ground state mode samples a large portion of the dipolar 
interaction. This softens the interaction, making it less 
repulsive than the centrifugai term at small p. When D is 
decreased, the softened interaction is pushed to larger p. 
li D is small enough the dipolar interaction entirely looses 
the competition with the attractive centrifugai term and 
there is no barrier. For coUisions in the first excited state 
(n = 1 and m = 0), Di = 0.95. This is important for 
identical for fermions, which will be able to collide in this 
channel without Wigner suppression. 



III. TUNNELING PROBABILITY 

We now investigate the tunneling behavior and some 
scattering properties of the quasi 2D system. To do this 
we assume that the scattering energy is smaller than 
the confinement energy and this means there is only one 
open channel (a unique combination of n and m). We 
solve Eq. © with the Johnson Log-derivative [20] and 
match to energy normalized Bessel functions. The cross 
section is (T = I sin^ (S) where S is the scattering phase 
shift. For pure threshold dipolar scattering the system 
behaves as a system with a 2D scattering length [21|| o f 
as/D = e2T+'"(2) _ 6.344, where 7 is Euler constantgl]. 
The m—O cross section diverges as energy goes to zero. 
Furthermore, in contrast to 3D the non-zero partial waves 
are Wi gner suppressed and go to zero as energy goes to 
zero 0,111. 

To explore the long range scattering of this system as 
a function of D and the scattering energy, we use the 
inner boundary condition to study the influence of the 
short range on the scattering. In future studies we envi- 
sion a full 3D solution combined with a trame transfor- 
mation ,23.] to capture the short range dynamics of the 




— Pure 2D s-wave 














(b) 



0.2 0.4 0.6 0.8 

Boundary Phasc/Ti 



1 10"- 



10" 10 

Sciittering Energy (E ) 



FIG. 3: (Color Online) (a) The scattering rate, 4sin^((5), as 
a function of the inner boundary phase for D=\ at different 
energies: 0.001 (solid red), 0.013 (black dashed), and 0.17 
(blue with circles). (b) The background scattering rate as a 
function of the scattering energy for D equal to 1 (red circles) , 
2 (open blue squares), and 5 (black triangle). The pure 2D 
scattering is shown as a solid blue curve. 



scattering. This would resemble the formalism used in 
quantum defect theory |24|. For this study, we simply 
place an inner wall behind the adiabatic barrier where 
the short range dynamics occur. A good example of the 
Information provided by varying inner boundary condi- 
tion is shown in Fig. [3] (a) . In this figure D — 1 and the 
energies are: 0.001 (solid black), 0.013 (red dashed), and 
0.17 (blue circles). The centrai feature is the shape reso- 
nance we want to analyze, and the other two resonances 
are narrow multi-channel Feshbach resonances. We use 
up to n = 8 to converge these scattering calculations. 

The scattering phase shift can be written as <5 w Sbg + 
tan^^ [rcot(^Si?, — '''■/2)] where Sbg is the back ground 
phase shift, F is resonance width, and 4>sR is the phase 
acquired inside the barrier. Sbg is the scattering phase 
shift acquired by the long range scattering and (t>sR is 
the impact of the short range dynamics on the scatter- 
ing. Here we directly alter 4'sr with the inner boundary 
condition. This leads to a direct method to estimate the 
width of the resonance produced by the dipolar barrier 
in a scattering process for a given D and energy. From 
curves like those in Fig. [3] (a) we are able to extract 
both Sbg and F as a function of 4>sr- We first look at the 
background scattering rate, which is shown in Fig. |3] (b) 
as a function of the scattering energy for D equal to 1 
(red circles), 2 (open blue squares), and 5 (black trian- 
gle). The important feature of this figure is that as D is 
increased, the cross section becomes more like the pure 
2D case (solid blue). Additionally as energy is lowered, 
the Systems becomes more like the pure 2D case. Re- 
stating the important point, for small D the scattering 
rate is smaller than the pure 2D scattering rate. This al- 
lows it to be used as an upper estimate of the scattering 
proprieties. 

We will now study the WKB tunnelling probability, 
and use this quantity to estimate of how likely the sys- 
tem is to access short range. Additionally the WKB tun- 
nelling probability can be used to estimate F 



251. R-eso- 



nances in this system will be connected to the threshold 
resonances which occur in dipolar scattering as the po- 
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larization is increased [2^ H^l • The confinement will 
modify the resonant character of the system, especiaUy 
when there is a barrier to the scattering process. As D 
is increased, the system gains another bound state and 
there will be a resonance. We can estimate F when the 
system is tunneling dominated. The WKB tunneling esti- 
mate will not be accurate when the system is scattering 
near the top of the barrier where much more complex 
scattering dynamic occur. 

The WKB tunnelling probability is: 



Ft = e ^"^,7 / dpK{p) 



(3) 



where k = y^2p{V^{p) — E), where is the low- 
est adiabatic potential including the Langer correction 
(m^ — j m^)- The integrai is over the classically for- 
bidden (c/) region for a given energy. In Fig. |4] (a) Pt is 
shown as a function energy for D equal to 1 (black line) , 
2 (red circles), 3 (blue triangles), 4 (maroon squares), 
and 5 (violet diamond). We also show the fitted widths 
from the scattering calculations for D =1 (black solid cir- 
cles), 2 (red open circles), and 3 (large blue triangle) as a 
function of energy. The important feature of this data is 
that the energy dependence of the resonance width is cor- 
rect; as the energy is decreased the width also decreases. 
One short coming of this analysis is that the tunneling 
width depends on the short range scattering. We try to 
avoid this impact by placing the hard wall near where the 
adiabatic potentials predict the barrier to end. However 
we only include the widths for _D — 1, 2, and 3 because 
for larger D the widths are very sensitive to the location 
of the inner boundary as channel couplings become very 
large. The resonances become more narrow the further 
back the wall is placed and requires many more channels 
to be converged. 

Here we provide an estimate of the WKB probability 
with a fit to the numerical data: 



In(PT) = aD"^'^" + hD^^/^° 
a = -5.17436 + 0.143167y- 



(4) 



0.00934332/-^ 



ò = 5.16437 + 1.61799y + 0.135513^^ 



where y = log{E / Ed). This fit is accurate to better than 
15% for D between 1 to 10 and energies between 10^® 
to 10^^. This was to capture the tunneling probability 
in the threshold regime and is much better for low en- 
ergy. Ref. [isj obtained an energy independent tunneling 

probability proportional to e~^ *^^ , the numerical Pt 
strongly deviates from this form as D becomes small and 
when the scattering energy is near the top of the barrier. 
This fit will be useful when estimating the nature of the 
scattering. 

In the semi-classical regime, E/Eo > 1, there are 
many partial waves scattering and the scattering rate ap- 
proaches 4%/ nDk 19]. The s-wave contribution is only a 
fraction of the scattering events. For non-zero m, the 
system has a repulsive centrifugai barrier in addition to 
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FIG. 4: (Color Online) (a) The WKB tunneling probability 
is shown as a function of energy for D equal to 1 (black line), 
2 (red circles), 3 (blue triangles), 4 (maroon squares), and 
5 (violet diamond). The fitted resonance widths from the 
scattering calculations are shown for D= 1 (black circles), 2 
(red open circles), 3 (blue open triangle). The green x's in 
each plot are the fitted Pt from Eq. (g]). (b) The WKB 
tunneling probability is shown as a function of energy for D 
equal to 1 (black line), 2 (red circles), 5 (violet diamond), 10 
(open blue circles), 20 (open black triangle), and 30 (open 
green squares). (c) Pt is shown as a function of D for the 
energies: 1 (black), 3.2 (red with circles), 10 (blue x), 32 
(green square), and 100 (violet diamond). The open green 



circles are the fit: 139.3(-^)^e 
100 and show very good agreement. 



=.2/5 



for E/Ed^W and 



5 




FIG. 5: (Color Online)The character of the quasi 2D scat- 
tering as a function of E/Ed vs. D. The excited threshold 
(solid blue) and the height of adiabatic barrier (solid red) are 
shown. The WKB tunneling probability contour for Pt ~ 
W'^ (dashed black), 10"^ (solid black), and 10"* (dash dot 
black) are plotted. The green x's in each plot are the fitted 
Pt from Eq. Q and the greens circles are the fit for the 
semi-classical energy regime from Eq. ([S]). 

the dipolar repulsion; this will greatly suppress the tun- 
nelling to short range. However to calculate the WKB 
tunneling probability, the short range interaction must be 
taken into account and for this reason it will be omitted 
here. 

As D is increased the impact of the s-wave reaching 
short range becomes much smaller. To illustrate the 
suppression of the tunnelling rates for m — 0, we have 
plotted Pt as a function of both energy and D in Fig. 
m (c). The WKB tunneling probability is shown as a 
function of energy for D equal to 1 (black line), 2 (red 
circles), 5 (violet diamond), 10 (open blue circles), 20 
(open black triangle), and 30 (open green squares). This 
figure demonstrates that as D is increased the tunneling 
rate will dramatically drop. The functional form of Pt 
is nearly CEe~^'^^^ ■ For D between 10 and 100 and 
energy between 10 and 100, we find: 

Pt - \m:i(ElEufe-^-^'>^"'" (5) 

This reproduces Pt to within about 30% over the stated 
parameter range. The fit is best when the energy is not 
near the top of the barrier. In Fig. |4] (c) this fit is shown 
as green circles for both E /Ed =10 and 100. Between 
the regimes when the energy is near the top of the bar- 
rier, the tunneling behavior is complex and will in fact 
involve much more than just tunneling dynamics to un- 
derstand. For this reason, we omit the parameter region 
from the fit. We have used up to n=8 and this is enough 
to converge Pt for the parameter range shown. 

Fig. [5] shows an estimate of effectiveness of the bar- 
rier and the character of quasi 2D scattering as a func- 
tion of E /Ed vs. D. The curves in this figure are: the 
first excited threshold (solid blue), the height of adia- 



batic barrier (red curve), the WKB tunneling probabil- 
ity contour for Pt = 10~^ (dashed black), 10~^ (solid 
black), and 10^** (dash dot black). Additionally in Fig. 
[S] the green x's and open circles are the contours from 
the fits, they agree well with the full data. The impor- 
tant feature of this figure is that it offers a quick estimate 
of whether or not a quasi 2D dipolar system will effec- 
tively produce a barrier to the short range. Once the 
parameters are below the dipolar barrier (red line), the 
system must tunnel to reach the short range interaction. 
When the parameters of the system are among or below 
the black contours the dipolar barrier will suppress the 
Systems ability to reach short range where inelastic col- 
lisions occur. More contours indicating higher levels of 
suppression were omitted because they offered little ad- 
ditional Information. In this case of highly suppressed 
tunneling systems, their behavior will be like the pure 
2D dipolar scattering, especially at threshold [T9I]. If the 
systems is near or above the dipolar barrier or if there 
is no dipolar barrier {D < De), the scattering dynamics 
will be rich and require the short range interaction to be 
included. For parameters above the solid blue line there 
will be many open confinement thresholds and the scat- 
tering dynamics will be very complex with many possible 
inelastic processes between confinement thresholds. 

In Fig. [S]the thin horizontal dashed line at E /Ed = 
1 denotes the transition between threshold and semi- 
classical scattering behavior. Above this line, many non- 
zero partial waves will contribute to the scattering and 
the cross section will behave as ka ex \/Dk (total scat- 
tering rate). When the scattering is in this regime, the 
s-wave scattering will be a fraction of total scattering 
events, and the larger number of non-zero partial wave 
scattering events might further suppress the inelastic 
losses. To know this for sure, one would have to include 
the short range interaction. 



IV. CONCLUSIONS 

In this paper we have studied quasi 2 dimensionai dipo- 
lar scattering. We included the effects of confinement by 
expanding this degree of freedom in harmonic oscillator 
States, and studied the properties of the resulting multi- 
channel system. We examined the adiabatic curves as a 
function of D, and found there is a dipolar barrier only 
when D > De = 0.34. We used the lowest adiabatic curve 
to obtain the WKB tunneling probability of this system. 
We found fits of Pt in the threshold and semi-classical 
scattering regime. When the system is tunneling domi- 
nated these will be good estimates of how likely particles 
are to make it to the short range. This might also be 
related to the width of resonances in the system. Fig. [5] 
offers a quick means to estimate whether an effective bar- 
rier will be produced by a molecular system in an optical 
lattice. 

The physical implications of this work are simple, try 

3/2^2^^1/2 

to maximize D = - — — . D must exceed 0.34 for 
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there to be a dipelar barrier, and should be much larger 
to significantly inhibit the particles from reaching the 
short range. If one desires to have a gas in the thresh- 
old regime, one must also keep E/Ed = ^ j^a^ less than 
1. This wih leads to an optimal value of d, which is set 
by the external electric field. To give a physical exam- 
ple, set E ~ lOOnK and lo ~ 50kHz. Then for LiCs, D 
can exceed 100. However the cost of this would be to 
make Ed minuscule, less than InK, and this would lead 
to semi-classical scattering. To maintain the thresholds 
scattering one would have set the field to d^E / B ~ 0.2 
where B is the rotational Constant of molecule and 
db is the bare dipole moment. At this field D ~ 2.7 and 
E / Ed ~ 0.2, which is in the threshold regime and pro- 
duces Pt ~ 2 • 10~^. An important physical example is 
RbK [11,0], for this system if the field is d^E/B ^ 5, the 
field will produce £> - 2, E /Ed ^ 0.1 and Pt ^i- 10"^. 
To further decrease Pt lowering the scattering energy 
might be more feasible than producing a tighter trap. 
This system will have threshold dipolar scattering and 
will help suppress the chemical reactions 0. There are 



several important future directions of this research. For 
example how does the scattering behave when E/huj > 2 
and there are many open confinement thresholds? How 
are the resonances in this quasi 2D systems related to the 
threshold resonances in 3D dipolar systems? To address 
this question a short range interaction must be included. 
This will also open up many avenues for further and more 
complete studies of these collisions. 
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